Geometric Microcanonical Thermodynamics for Systems with First Integrals 
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In the general case of a many-body Hamiltonian system, described by an autonomous Hamiltonian 
H , and with K > independent conserved quantities, we derive the microcanonical thermodynam- 
ics. By a simple approach, based on the differential geometry, we derive the microcanonical entropy 
and the derivatives of the entropy with respect to the conserved quantities. In such a way, we show 
that all the thermodynamical quantities, as the temperature, the chemical potential or the specific 
heat, are measured as a microcanonical average of the appropriate microscopic dynamical functions 
that we have explicitly derived. Our method applies also in the case of non-separable Hamiltonians, 
where the usual definition of kinetic temperature, derived by the virial theorem, does not apply. 
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In spite traditionally, the microcanonical ensemble has 
played a minor role respect to the canonical and gran- 
canonical approach, in reason of the more complicated 
formulation of the statistical formulas, nowadays, thanks 
to the performance provided by the current computers, 
the microcanonical simulation is probably the most used 
tool in physics. Molecular dynamics is applied for investi- 
gations that cover the fields of the dynamics of nonlinear 
semi-classical systems, the living matter issues (protein 
folding, DNA- helix), the dynamics of nanosystems and 
so on. The ergodicity in fact, whenever holds true, al- 
lows to measure dynamic observables as temporal aver- 
ages of appropriate functionals along the dynamical evo- 
lution, instead of perform integrals (often unmanageable) 
on the phase space. In this contest have been written 
very famous papers [l|, that date back to the early ap- 
plications of molecular-dynamics in physics. However, 
their application is restricted to systems described by a 
separable Hamiltonian, i.e. given by a standard term of 
kinetic energy plus an interaction potential term. In a 
series of more recent papers [H-Q: Rugh has developed 
a microthermodynamic formalism to measure thermody- 
namical functions in the microcanonical ensemble. The 
Rugh's formalism applies to non-separable Hamiltonian 
systems, and to the case of Hamiltonians that depend on 
parameters. 

In the present paper, we illustrate a simple geo- 
metric approach to measure thermodynamic observables 
within the microcanonical ensemble, in systems described 
also by nonstandard Hamitonians Q and with an ar- 
bitrary number of first integrals. We consider a clas- 
sical many-particle system described by a Hamiltonian 
Hix 1 , . . . , x N+K ) which, in addition to the total en- 
ergy H, has K > independent conserved quantities 



such that {H,Vj} = for j 



1, . . . , K. All the first integrals are assumed to be in invo- 
lution. We derive in detail the expressions for the micro- 
canonical entropy (i.e. the expression of the microcanon- 
ical invariant measure ), and we give the general formulas 
that allow to measure all the thermodynamic function as 
the temperature, the chemical potential, the specific heat, 



the pressure, and so on. More precisely, we show that 
entropy and the other thermodynamic functions are de- 
rived by multidimensional integrals performed over the 
sub-manifolds given by the intersection of the constant 
energy hyper-surfaces H(x 1 , . . . ,x N+K ) = vq with those 



defined by Vj (x 1 , . 



for j 



1, 



,K. In 



particular, we show that all order derivatives of the en- 
tropy with respect to the Vj for j = 0, . . . , K, from which 
the microcanonical observables depend on, can be calcu- 
lated as time averages along the time evolution of appro- 
priate functions, whenever the hypothesis of ergodicity 
holds true. In order to simplify the notations, let us intro- 
duce the following definitions: with x := (x 1 , . . . , x N+K ) 
we indicate a generic point of the system phase-space, 
with Vb(x) := H(x) the "first" of the conserved quan- 
tities, and with Aj. = {x e R N+K \ Vj(x) = Vj} for 
j = 0, . . . , K the hyper-surfaces associated to the K + 1 
conserved independent quantities. We shall assume the 
level sets of Vj to be non-singular hyper-surfaces. Even 
if the hyper-surfaces A J V , in general, constitute a singular 
foliation, since for some values of Vj an hyper-surface 
is not a differential manifold, for generic values of Vj 
this is not an issue. Furthermore, we will indicate with 
V = (Vq,...,Vk) a vector whose components are the 
conserved quantities, and with v = (vq, . . . , vk) the vec- 
tor of the corresponding values, for a generic configura- 
tion. In addition, we will indicate with OTjr(v) C K N+K 
the set given by the intersection of the hyper-surfaces 
Aj . , for j = 0, . . . , K, that is DJt K (v) = nf =0 K 3 ■ We 
shall give the proof of statements reported above, in three 
steps, according to the choice K — 0, 1 and for an arbi- 
trary value of K. 

Case K = 0.— Let us begin by summarizing the re- 
sults of Ref. Let #(x) be a classical Hamiltonian 
describing an autonomous many-body system whose co- 
ordinates and canonical momenta are represented as N- 
componcnt vectors x G M. N . Let us suppose that this 
system has just one conserved quantity, that is the total 
energy H. The microcanonical description of equilibrium 
thermodynamics for this system, is given in terms of the 
microcanonical entropy S(vq). Among the equivalent ex- 
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pressions allowed for S(vq) one can consider the surface 
entropy [6] S(v ) = \nfd N xS(V (x) - v ), in @, i] it is 
shown the following identity 



where F„ 



S(v ) = In 



A" 
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This expression has a precise geometrical interpretation. 



In fact, m N (Vb) 
is the metric induced 
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A?,., where nS = 9,-Vb 
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dxj ■ ■ • dxr 
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on the hyper-surface 

,, , ru/ 1| v roll are the components of 

the unitary vector orthogonal to the hyper-surface, and 
l/||VVb|| is the microcanonical measure. The symbol 
dxj, means that dxj has been lifted from the formula. 
The microcanonical thermodynamics is obtained by cal- 
culating the derivatives of S(vq). This fact allows a 
geometric interpretation of all thermodynamic quantity. 
E.g. the inverse temperature is given by the definition 
1/T(vq) = dS(vo)/dvo- This quantity can be rewrit- 
ten in a geometric form by using the Federer-Laurence 
derivation formula 



m N (Vb) V = / m N (Vb) A k (Vb, r/>) (1) 



gk 

where, by setting n° = VVb/|| VVb||, it is 

A(V ,i>) = 1/|| W ||V (n» . (2) 

In fact, the inverse microcanonical temperature results 
(see [|[ for details) 

1 X A om*(Vb)/||W ||V(n /||W||) 

"0 



T(vo) 



Jao m"(Vo)/||VVb|| 



that is, the microcanonical average of V (n° /\\ Wb||) on 
the energy hyper-surface A° o . 

Case K = 1.— Let us consider a system described by 
the Hamiltonian H (x 1 , . . . , x N+1 ) which, has in addition 
to the total energy iJ, one independent conserved quan- 
tity V 1 (x 1 , . . . , = vi such that {H, Vi} = 0. Con- 
sistently with the previous case, the microcanonical en- 
tropy is given by S(v) = lnJ'd Ar+1 x<5(V(x) — v), where 
<5(V(x)-v) = Uj=o,i S ( v j( x )~ v j)- The second delta 
function has been added in order to take in account of 
the second first-integral of motion. In [12| it has been 
shown that 



S(v) = In f 



m ™(V ,Vi) 
IHiII 



where 9Jti is the intersection set of the hyper-surfaces A° o 
and A^ . The measure m N+1 (Vb, Vi) is that one induced 
by R^ 1 on Tti, that is 



N+l 



m 



,+1 (Vo,^i) F flu dx 1 ...dx tl ...dx u ...dx 

(rz 1 ) 



N+l 



-e° e* ) depends on the unitary vectors 
of a basis {e^e 1 } are derived by the normalized gradi- 
ent vectors n? — W}/||VV}|| for j = 0, 1, by means of a 
Gram-Schmidt orthonormalization process, starting from 
the vector e 1 = nr. Thus, e° = [n° — (n° • ft 1 )?! 1 ]/!! — 
(n° ■ n 1 ) 2 ] 1 / 2 . Furthermore III = Wq A VVi, where 
VVj — ^Z^d^Vjdx^ for j = 0, 1. In order to derive 
the temperature in the microcanonical ensemble, accord- 
ing to the definition T(v) = (dS(vo, vi)/dv ) ~ , we shall 
use the following generalization of the Federer-Laurence 
derivation formula ([l])-©, whose proof is given below. 
The generalized derivation formula results |13| 

|^/m- 1 (V ,Vi)^=/ m^(V ,V 1 )A h ;(V ,V 1 ,i;) (3) 
where 



Ai(Vb,Vi,V) 



^[V^VO-^-^-V) (e% 

(4) 

is meant as a function of Vb and Vi also through the 
dependence from these latter of the unitary vectors e° 
and e . Here, and in the following we mean e 1 • (e 1 ■ 
V)(e°) = E^e}4d fe (e°). 

By the calculation reported above, we can easily de- 
rive dS(E, v\)/dvi that, when Vi is the total number 
of atoms, gives the chemical potential via the definition 
fx/T = — dS(E, vi)/dv\. Indeed, by exchanging Vb with 
Vi and vq with v\, it is easy derive from Eq. ([3]) the 
following formula 



m N+1 (Vb,W 



m N+1 (Vb,Vi)^ (Vb,Vi,^) 



strti 



where 



Obviously 



B 1 (Fo,F 1 ^)=A 1 (Vi,Fo,V)- 



(5) 



gk+q 



m N+1 (V ,Vi) V 



m^{V ,V{) Bf (Vb.Vi.Af (Vb,Vi») • 



(6) 



Proof of the generalization of the Federer-Laurence 
derivation formula.— hi this section we shall give the 
proof of the generalization of the Federer-Laurence the- 
orem to varieties of co-dimension two. As first step we 
shall calculate d(f m vpm N+1 (Vo, V\))/dvo and, by itera- 
tion we will achieve our aim. In this derivation we make 
use of the Stokes theorem. Let U be a subset of A*, 
such that dU = VJtx{v + Av ,v)umi(v ,v). The Stokes 
theorem states 



u = / dui . 

du Ju 



(7) 
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Let us choose ui = -0m N+1 (Vq, Vi), thus we obtain 



JV+l 



duj = ^2 {—l) u ®fi {'ipFfii') dxi . . . dx v . . .dx N+ i 

1 

and, by using the definition of F^ Vl it follows 



dcj = / [V • (ipe ) - ^e 1 • (e 1 • V) (e )] m N+1 (Vi) . 
u Ju 



Now it is easy to verify that 

m N + i (vi) = e° a dx a A m^ 1 (V , Vi) , 

a 

therefore, by choosing vectors orthogonal to WIi(vq,v) 
and 9JXi(vq + Avq,v) (but tangent to A*) with opposite 
orientations, we get from the Stokes theorem 

/ Vm K+1 (F ,^i) - / ^m^ 1 (V Q ,V 1 ) = 

Jwti(v +Ava,v) Jwii(v ,v) 
Av p 

ds / m K+1 (Vb, Vi) x 

o Jmi(v +s,v) 

[V • tye°) - ^ ■ (e 1 • V) (e )] , 

that, in the limit Avq —> gives Eq. © for k = 1, 
by iteration the generalization of the Laurence-Federer 
formula is proved. 

General case.— In this section we will derive the micro- 
canonical entropy, and its derivatives with respect to the 
conserved quantities, in the case of systems with K > 1 
first integrals in addition to the energy. 

Let a system be described by the Hamiltonian H(x) 
that depends on x = (x 1 , . . . , x N+K ). Let us suppose 
the system to have, in addition to the total energy H, 
K independent conserved quantity V,(x) = Vj such that 
{H, Vj} = for j = 1, . . . , K. The microcanonical en- 
tropy is given by S(v) — \nJd N+K xS(V(x) — v), where 

5(V(x)-v) = Yl?=a 5 ( v j(x-)- v j)- In this case 



S(v) 



In 



an. 



m K + K (v) 



where ^XSIk is the intersection set of the hyper-surfaces 
Aj. for j = 0,...,K, V = (V ,...,V K ) and v = 
(vq,...,vk)- The measure m K+K (V) is that one in- 
duced by R N+K on 9JIk, and it is defined as follows. 
Let n? — VVj/||VVj|| be the unitary vectors orthogo- 
nal to the hyper-surfaces Aj., for j = 0, . . . , K. By 
the Gram-Schmidt orthonormalization process, start- 
ing from the set of K + 1 independent unitary vectors 
{n°, . . . , n K } we can obtain the following orthonormal- 
ized basis {e°, . . . , e }, where e K — n K . Therefore, the 
measure m N+K (V) results 



where e- 7 = ^ a e^di: <> , for j = 0,...,K, and * is the 
(Hodge) star operator. The microcanonical measure is 
the norm of 



U K = W A Wi A • • • A VV K , 



(8) 



where Wj = ^ d^Vjdx^, for j = 0, . . . , K. 

For convenience let us introduce the following nota- 
tion: (v) = (v\, . . . ,Vk) and (V) = (Vi, . . . , Vr-). Let 
us calculate the derivative of entropy with respect to 
«o, the derivates with respect to the other variables 
are obtained easily by cyclic permutation of the indices. 
Also here we resort to the Stokes theorem (J7]) where, in 
this case, U is a sub-set of the intersection set of the 
hyper-surfaces A? . for j = 1, . . . , A , in such a way that 
dU = 9Rk{vo + Au , v) U %JIk{vo, v). Furthermore, we 
choose to = ^ntl N+K (V), thus the following facts hold true 

m N+K (V) = e *dx a A m N - K (V) , 

OL 

where m N+K (V) = *(e A ' A • • • A e 1 ) 



/ tpm N+K (V)- ^m N+K (V) = 

Jm K (v +Av ,v) J<m K (v ,v) 

ds / m N+K (V) x 

ii Jmi(v a +s,v) 



^(dw.^A-Ae 1 )) 



VVo-e 

where (, ) is the inner product. Thus we have 
d< 



m N + K (V)V=/ m N + K (V) A^(V,V) (9) 
UV Jvji k Jsut K 



where 

A K (V,iP) 



1 



Wo ■ e' 



^(du,*(e K A-'-Ae 1 )) . (10) 



m N+K (V) = *(e A A 



Ae u 



Application to the discrete nonlinear Schrodinger 
equation.— With reference to the standard canonical co- 
ordinates, the Hamiltonian of the discrete nonlinear 
Schroringer equation (DNSE) writes 

H a (p, q) = A a J2{p 2 j + q 2 jf ~ r a ^2(pjPj +1 + qjq 3+l ) , 



where A Q and r Q measure the magnitude of the repulsive 
on-site interaction and of the hopping intensity, respec- 
tively. Here, the index j numbers the sites from 1 to 
M and periodic boundary conditions are assumed. This 
Hamiltonian describes the dynamics of a system of a 
Bose-Einstein condensate (BEC) of repulsive-atoms, in 
an optical lattice in the superfluid regime. The model 
possesses two conserved quantities: the energy Vo = H a 
and the number of particles Vi = YljiQj Thus, 
it corresponds to the case K — 1 of the present paper. 
Under the hypothesis that the ergodicity holds true, for 
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this system the temperature, as a function of the total 
energy E and number of atoms N, can be measured as 
a temporal average along the dynamical evolution. By 
using Eq. (U]) we get 



T(E,N) s 



lim — 

oo s 



where 



V 



III 



ds'Mx(s'))} , 



N) (e°) 

liriiii 



is written in terms of e 1 = Wi/||VVi|| and 

e = [n -(n°.e 1 )e 1 ]/[l-(» -e 1 ) a ] 1/a » 
where n° = VH a /\\S/H a \\ . Further more it results 

1/2 



in. 



v 1 ■' 'nu., dVx OJi a <)v\ N 



E 



dxp dx u dx v dx^ 



with and x^-i — q^, x 2fJ , = Pn are the M = (N + l)/2 
lattice coordinates. By exchanging H a with Vi, and vice 
versa, in these equations we can measure the chemical 
potential. In Ref. 1J| , these formulas have been checked 



and used to explore the thermodynamics of the DNSE 
where it has been observed the existence of negative tem- 
perature states. Analogously, in Ref. [l5| the micro- 
canonical approach proposed by Rugh, has been applied 
in the contest of Bose-Einstein condensates, studying the 
projected Gross-Pitaevskii equation. 

Application to the BEC mixtures.— The dynamics of a 
mixture of two bosonic species, condensed in an optical 
lattice and in the superfluid regime is described by two 
coupled DNSEs and the Hamiltonian writes 



H(p,q;p,q) = H 1 (p,q) + H 2 (p,q)+f3j2(p 2 j+qj)(P. 



* ' 9?) 



where /3 takes into account the interaction between the 
two species, and by assuming periodic boundary con- 
ditions the index j runs from 1 to M. In this case 
there are three conserved quantities the total energy 



Vq = H. and the total number of atoms of each species: 
Vi = Y,j(q] +J#/2 and V 2 = Y^(gf +P-)A Thus, we 
have to use the formulas of the case K = 2 and, from Eq. 
(|10p the inverse temperature is measured as the micro- 
canonical average of the function 



in 



WH-e° 

where S is the Kronecker's delta and 

N+K 2 
p=l i,j,k=Q 

Here, e is the completly antisymmetric tensor, the uni- 
tary vectors are = VVj/||VV^|| for j = 1,2, and 
e° = 77/||f?||, where 

r) = n° - (n° • e l )e l - (n° • e 2 )e 2 , 

and n° = VfT/||Vif||. Finally, from Eq. © we have 



in 2 || = 



t J >i<t*j<t*k 



dH dVi dV 2 
dx^ dx H dx^ k 



1/2 



with 2^-1 = q^, x 2fJ . = p^, for fj, = 1, . . . , (N + 2)/4, 
and x 2m -i = g^-M, ^2 M = P^-Af , for /i = (N + 2)/4 + 
1, . . . , (iV + 2) /2. In this case the number of lattice sites is 
M = (N+2) /4. These equations are presently used to in- 
vestigate the thermodynamics of two-component bosons 
mixtures in Ref. (l6| . 

Final remarks.— In this paper we have proposed a geo- 
metric approach to measure thermodynamic observables 
within the microcanonical ensemble in classical Hamilto- 
nian systems with K > independent conserved quanti- 
ties. The method that we have shown is applicable also 
to systems with nonstandard Hamiltonians and with an 
arbitrary number of conserved quantities. As an exam- 
ple we have derived the formula of the temperature for 
a model that describes the superfluid dynamics of two- 
components bosons on an optical lattice. 
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